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Abstract 

A search for neutral Higgs bosons has been performed with the OPAL detector at LEP, 
using approximately 170 pb~^ of e+e" collision data collected at ~ 189 GeV. Searches 
have been performed for the Standard Model (SM) process e~^e~—>-}i^Z^ and the MSSM 
processes e"'"e~— i>h°Z°, A°h°. The searches are sensitive to the bb and r+r~ decay modes 
of the Higgs bosons, and also to the MSSM decay mode h°— ^A'^A'^. OPAL search results 
at lower centre-of-mass energies have been incorporated in the limits we set, which are 
valid at the 95% confidence level. For the SM Higgs boson, we obtain a lower mass bound 
of 91.0 GeV. In the MSSM, our limits are mh> 74.8 GeV and mA> 76.5 GeV, assuming 
tan/?>l, that the mixing of the scalar top quarks is either zero or maximal, and that the 
soft SUSY-breaking masses are 1 TeV. For the case of zero scalar top mixing, we exclude 
values of tan (3 between 0.72 and 2.19. 
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1 Introduction 



The OPAL detector at LEP has collected more than 185 pb~ of e"'"e~ collision data 
at centre-of-mass energies in the vicinity of 189 GeV. These data are used to search for 
neutral Higgs bosons [|1| within the framework of the Standard Model (SM) 0, and within 
the Minimal Supersymmetric extension of the Standard Model (MSSM) 0. Searches have 
been performed for the "Higgs-strahlung" process, e'^e~^h.^Z^^h.^f{, where h° is either 
the SM Higgs boson Hgp^ or the lightest neutral CP-even Higgs boson in the MSSM, and 
ff is a fermion-antifermion pair from Z° decay. For the h^uu (h°e^e~) final state, the 
contribution from the WW (ZZ) fusion process, which is only a small part of the total 
production except near the kinematic threshold of the e^e~— >h°Z° process, is also taken 
into account. We have also searched for the MSSM process e+e~— i>A°h°, where A° is the 
CP-odd Higgs boson. 

In this paper, only the dominant decays of the neutral Higgs bosons (h° and A°) into 
bb pairs and r^r~ pairs are considered. We also consider the MSSM decay h°— i>A°A° 
with A°^bb. Searches for Higgs boson decays into SUSY particles are not presented in 
this paper. OPAL searches at centre-of-mass energies up to 183 GeV have been reported 
in 0, 1^, 1^ for the neutral Higgs bosons of the SM and the MSSM. For the SM Higgs 
boson a lower mass bound of mh > 88.3 GeV was obtained at the 95% confidence level 
from the lower-energy searches. The searches presented here are similar in procedure to 
our previous searches at ^/s ^183 GeV [Q; in this paper only significant changes from our 
previous analyses are described in detail. In particular, we have significantly improved 
our b-tagging algorithm; the new algorithm is described in detail in Section ^ 

We have optimised the selection requirements in each of our searches separately to 
maximise the sensitivity to Higgs boson production. In each case, the mass limit we would 
expect to set in the absence of a Higgs signal, computed using Monte Carlo simulations 
of signal and background, has been maximised. When the distribution of a variable in 
the signal Monte Carlo is needed as an input to a selection technique, we use a mixture 
of Higgs masses near the expected limit in order to optimise the analyses over the range 
of Higgs masses under study. 

Recent searches for neutral Higgs bosons performed by the other LEP collaborations 
are listed in The combined mass limit by the four LEP collaborations for the SM 
Higgs boson and a combination of the MSSM Higgs boson searches using data taken at 
^/s < 183 GeV are reported in [Q. 

2 OPAL Detector, Data Sets and Monte Carlo Sam- 
ples 

The present analysis is based on data collected with the OPAL detector during 1998 at 
a luminosity-weighted average centre-of-mass energy of 188.6 GeV. The searches presented 
here use subsets of the data sample for which the necessary detector components were 
fully operational. Approximately 170 pb^^ were analysed, varying ±2% from channel to 
channel, depending on the subdetectors required. 

The OPAL detector has nearly complete solid angle coverage. The central tracking 
detector consists of a high-resolution silicon microstrip vertex detector which im- 
mediately surrounds the beam pipe. Its coverage in polar angle[] is |cos0| < 0.9. The 

^ OPAL uses a right-handed coordinate system where the +z direction is along the electron beam 



silicon microvertex detector is followed by a high-precision vertex drift chamber, a large- 
volume jet chamber, and z-chambers to measure the z coordinate of tracks, all in a 
uniform 0.435 T axial magnetic field. The lead-glass electromagnetic calorimeter with 
a presampler is located outside the magnet coil. It provides, in combination with the 
forward calorimeters, the forward scintillating tile counter (the "MIP plug") [0, and the 
silicon-tungsten luminometer [O , geometrical acceptance down to 25 mrad from the beam 
direction. The silicon-tungsten luminometer serves to measure the integrated luminosity 
using small-angle Bhabha scattering events |]T3 . The magnet return yoke is instrumented 



with streamer tubes and thin gap chambers for hadron calorimetry; it is surrounded by 
several layers of muon chambers. 

Events are reconstructed from charged-particle tracks and energy deposits ("clusters") 
in the electromagnetic and hadron calorimeters. The tracks and clusters must pass a set 
of quality requirements similar to those used in In calculating the total visible 

energies and momenta, i?vis and Pvis, of events and individual jets, "energy-flow objects" 
are formed from the tracks and calorimeter clusters [|T5|. A track and its associated clusters 



are grouped into a single energy-flow object after correcting the clusters for the estimated 
energy deposited by the particle leaving the track. Unassociated clusters are left as energy- 
flow objects. This procedure prevents double-counting of the energies of particles which 
leave signals in more than one subdetector. 

A variety of Monte Carlo samples has been generated in order to estimate the detection 
efficiencies for Higgs boson production and SM background processes. Higgs production 
is modelled with the HZHA generator |jT6| for a wide range of Higgs masses. The size of 
these samples varies from 500 to 10,000 events. The background processes are simulated, 
with typically more than 50 times the statistics of the collected data, by the following 
event generators: PYTHIA ^ ((Z/7)*^qq(7)), grc4f pf and EXCALIBUR p| (four- 
fermion processes (4f)), BHWIDE (e+e-(7)), KORALZ (gg (/i+/i-(7) and r+r-(7)), 

Ml (hadronic and leptonic two-photon 



and PHOJET HERWIG |2|, and Vermaseren 
processes (77)). The hadronisation process is simulated with JETSET with parame- 
ters described in ||2^. The cluster fragmentation model in HERWIG is used to study the 
uncertainties due to fragmentation and hadronisation. For each Monte Carlo sample, the 
detector response to the generated particles is simulated in full detail [EH]. 



3 Improved b- Tagging 

The predominant decay mode of the SM Higgs boson in the mass range under study is 
expected to be to pairs of b quarks; the same is true of the neutral Higgs bosons of the 
MSSM in large areas of the model's parameter space. Therefore, efficient and pure tagging 
of b quarks is an essential technique in searches for these particles. We have developed a 
tagging method that uses the three nearly independent techniques of lifetime, high-pt 
lepton and kinematic tagging to identify jets containing b hadron decays and to minimise 
the contamination from jets which do not. Artificial Neural Networks (ANN's) have been 
introduced to combine optimally several lifetime-sensitive tagging variables and also to 
combine kinematic variables in the jet-kinematics part of the tag. For each jet, the outputs 
of the lifetime ANN, the kinematic ANN, and the lepton tag are combined into a likelihood 
variable B which discriminates b-flavoured jets from c-flavoured and uds-flavoured jets. 



and where +x points to the centre of the LEP ring. The polar angle Q is defined with respect to the +^ 
direction and the azimuthal angle (/> with respect to the -\-x direction. 



While the basic scheme of the tagging procedure has been kept unchanged with respect 
to our previous procedure 0, a significant improvement in the performance has been 
achieved for this analysis by introducing input variables with greater sensitivity. 

In order to improve the performance of the lifetime-sensitive part of the tag, we have 
developed a new algorithm to identify displaced vertices. An ANN applied to individual 
tracks (track- ANN) has been trained, using Z'^-pole Monte Carlo samples, to discriminate 
between tracks from the primary vertex and those from secondary vertices. It uses, among 
other inputs, the impact parameter of the track with respect to the primary vertex and 
the transverse momentum with respect to the jet axis. The tracks belonging to a given 
jet are ranked in descending order of their probability of having come from a secondary 
vertex, according to the track- ANN output. Using the first six tracks (or all tracks in the 
jet if there are fewer than six), a "seed" vertex is formed using the technique described 
p^ . This technique first forms a vertex from the input tracks, and then removes the 



m ^ ^ 

track which contributes the most to the vertex x^- This procedure is repeated until no 
track contributes more than 5 to the x^, and the seed vertex is the common fit origin of 
the tracks. The remaining tracks in the jet are then tested to see if they may be added 
to the seed vertex, based on their contributions to the vertex x^- 

In addition to identifying displaced vertices, we introduce a new method for combining 
track impact parameters in order to gain b-tagging efficiency in events where the secondary 
vertices are less distinct. In this method, the impact parameter significances S^'^ and 
S""^, in the r0 and rz projections, respectively, are formed by dividing the track impact 
parameters by their estimated errors. The distributions of S^'^ and S"'^ for each quark 
fiavour from Z°-pole Monte Carlo are used as the probability density functions (PDF's) 
f^"^ and /^^ (g=uds, c and b). The combined estimator Fg for each quark fiavour q is 
computed by multiplying the f^'^ and /^^ for all tracks passing the quality requirements. 
The final estimator Cjp is obtained as the ratio of Fj, and the sum of Fuds, Fc and Fi,. 

The following four variables are used as inputs to the lifetime ANN. They replace the 
five variables used in our 183 GeV analysis 

• The vertex significance likelihood {Cy): The likelihood for the vertex significance 
is computed analogously to the Cip above, using the decay length significance of 
the secondary vertex rather than the impact parameter significance of the tracks. 
Because the decay length significance of a secondary vertex depends strongly on the 
number of tracks in it, the PDF's for each fiavour are computed separately for each 
value of the charged multiplicity of the secondary vertex. 

• The reduced significance likelihood (JZy)'- To reduce sensitivity to single mismea- 
sured tracks, the track with the largest impact parameter significance with respect 
to the primary vertex is removed from the secondary vertex candidate and the re- 
maining tracks are used to recompute the likelihood £v If the original vertex has 
only two tracks, the function is calculated from the impact parameter significance 
of the remaining track. 

• The combined impact parameter likelihood {Cip) described above. 

• The reduced impact parameter likelihood iTZip): The track having the largest impact 
parameter significance has been removed in the calculation of £ip. 



In the jet-kinematics part of the tag, our previous procedure used only the boosted 
sphericity of the jet as an input to the final likelihood. This variable has been replaced 



with three new variables, which are combined with a separate ANN. These three inputs 
are 1) the number of energy-flow objects around the central part of the jet, 2) the angle 
between the jet axis and its boosted sphericity axis, and 3) the C-parameter [^] for the 
jet boosted back to its rest frame. The high-p^ lepton tag has not been changed. 

The outputs from the lifetime ANN, the jet-kinematics ANN and the high-pf lepton 
tag are combined with an unbinned likelihood calculation as described in , and the final 
output B is computed for each jet. Figure 0(a) shows the distribution of B for calibration 
data taken at ^/s = rrizo in 1998; it also shows the improved performance of the new b- 
tagging algorithm, compared with the previous version . For the same efficiency, the new 
b-tag has about two-thirds of the background acceptance of the previous tag. The tagging 
efficiency for b-flavoured jets has been verified to an accuracy of 1% with a double-tagging 
technique using the Z°-calibration data collected with the same detector configuration and 
operating conditions as the high-energy data; a comparison of the results of this double- 
tagging test between the data and Monte Carlo is shown in Figure 0(b). Using the 
high energy data, the expectation from the SM Monte Carlo agrees very well, within a 
relative statistical uncertainty of 5%, for samples enriched in (Z/7)*— >■ qq processes with 
and without hard initial-state photon radiation. For the lighter fiavours, the efficiency 
has been examined by vetoing b flavour in the opposite hemisphere, and the resulting 
estimate of the background is found to be described by Monte Carlo within an accuracy of 
5-10%. This constitutes one component of the systematic uncertainty in the background 
estimates of the search channels presented here. The efficiency for tagging lighter flavours 
has also been checked by computing B for a high-purity sample of light-flavour jets in 
W"'"W~ qqiu decays; there is no evidence of mismodelling of the light-flavour tag rate 
within the statistical precision of the test. 



4 Searches for e+e — ^h^Z° 

Throughout this section, h° denotes both the Standard Model Higgs boson Hg^^ and the 
lightest CP-even Higgs boson in the MSSM. We have searched for the process e"'"e~— s>h'^Z° 
in the following flnal states: h'^Z'^^bbqq (the four-jet channel), h°Z°— s>bbz/i/ (the missing- 
energy channel), h°Z°^bbr+r^ and r+r^qq (the tau channels), h^Z*^— >-bbe"'"e~ and bb/i"*"/i~ 
(the electron and muon channels). The selections for all channels are similar to those 
described in 0; their main features and differences from our previous publication are de- 
scribed in Sections [4. 1| through |4.4|. These selections are also sensitive to the h°Z°— >A°A°Z° 



processes of the MSSM, and the procedure for incorporating this additional information 
into our limits is described in Section [4.5| . 

4.1 The Four- Jet Channel 

The selection for the four-jet channel has been modifled with respect to our previous 
publication in order to improve the sensitivity to Higgs signals at higher masses. In 
addition to using the improved b-tagging algorithm described above, we have improved 
the grouping of flnal-state particles into four jets and the assignment of pairs of these jets 
to the Z" and h° bosons. 

The correct assignment of particles to jets plays an essential role in reducing one of 
the main backgrounds, W^'^W^^qqqq, and also in accurately reconstructing the Higgs 



boson mass. The Durham algorithm [29| is used to flnd four jets in each event, i.e., the 



resolution parameter ycnt is chosen to be between ^34 and ^45, where ^34 is the transition 
point from three to four jets, and 7/45 is the transition point from four to five jets. These 
jets are used as reference jets in the following procedure. Each particle is reassociated to 
the jet having the smallest "distance" in the EO scheme |Q of the JADE algorithm 



to the particle. This "distance" is defined to be Ej^^ ■ Eparticie ■ (1 — cos^,), where Ej^^ 
is the energy of the i*^ reference jet {i=l to 4), Eparticie is the energy of the particle 
(charged track and/or calorimeter cluster, corrected for double-counting) and 6i is the 
angle between the i^^ reference jet and the particle. This reassociation procedure reduces 
the fraction of wrong association of particles to jets, which improves by about 10% the 
di-jet mass resolution for the W^W~ background as well as for Monte Carlo Higgs signals, 
before kinematic fitting. The rearranged jets are used in cuts 5 and 6 of the preselection 
and the final likelihood selection described below. This procedure is also used for the 
four-jet h°A° channel described in Section O. 



The preselection is exactly the same as that used in ||^; it is designed to retain only 
four-jet-like events. The requirements are: 

1. Events must satisfy the hadronic final state requirement of Reference 



2. The effective centre-of-mass energy y s' [32] must be at least 150 GeV. 



3. The value of 1/34 in the Durham algorithm must exceed 0.003. 



4. The event shape C-parameter |28| must be larger than 0.25. 



5. Each of the four jets must have at least two tracks. 

6. The probabilities must be larger than 10~^ both for a four-constraint (4C) kine- 
matic fit, which requires energy and momentum conservation, and a five-constraint 
(5C) kinematic fit, additionally constraining the invariant mass of one pair of jets 
to rriio The algorithm used to select the jet pair which is constrained to the Z° 
mass is described below. 

There are six possible ways to assign four jets in pairs to the Z° and the h°. An 
algorithm based on a likelihood technique has been developed in order to select the most 
Higgs-like combination. This algorithm makes use of b-tagging variables and the proba- 
bility that the 5C kinematic fit would have had the observed or larger. For each of 
the six jet-assignment combinations, we form two variables from the b-tagging output as 
follows: (1) Bi ■ B2 for jets assigned to the Higgs boson, and (2) (1 — B3) ■ (1 — B4) for jets 
assigned to the Z°. These two variables, together with log(p5c) for the jet combination, 
are used as inputs to the likelihood calculation, which is designed to distinguish between 
the three cases of correct jet assignment, swapped pairings in which the two jets from the 
h° are wrongly assigned to the Z°, and other combinations, for which one jet from the 
h'' is assigned to the Z°. For each of the possible assignments of pairs of jets to the Z° 
and the h° with p^c > 10~^ three likelihoods, c^wapped ^others^ calculated 

using the PDF's of the input variables in signal Monte Carlo events passing the prese- 
lection. The relative likelihood £ = H^^^''^^ j (^c^owect ^ £s«,apped ^ £other.^ ^j^^^ formed 

for each combination and the combination yielding the largest value of L is chosen. The 
reconstructed Higgs mass is taken to be the mass obtained from the 5C kinematic fit for 
this combination. 



With this algorithm, the rate of correct jet pairings for selected signal events with 
TJih = 95 GeV is 71%, compared to the rate of 36% obtained with the method used in 
our 183 GeV analysis, which uses only p^c for the assignment. This improvement in 
the correct jet-pairing reduces the fraction of combinations for which the reconstructed 
h° and Z° masses have been swapped. The fraction of signal Monte Carlo events with 
TUh = 95 GeV for which the reconstructed Higgs mass is in the Gaussian core centered 
on 95 GeV increases by 10% when this improved jet pairing is used, as compared to 
the previous pairing. The correct jet pairing provides better assignment of the b-tag 
variables to the Higgs candidate and improves the separation of the h^Z" signal from Z°Z*^ 
background. 

After the preselection, the eight variables described in |^ are combined using a likeli- 
hood method P] . The distributions of four of these variables are shown in Figure |^ for the 
OPAL data and corresponding SM background simulations. The likelihood distribution 
C^^ and the distribution of the reconstructed Higgs boson candidate mass are shown in 
Figure ^ Higgs boson candidate events are required to have > 0.96. The numbers 
of observed and expected events after each selection step are given in Table |l|. 

After the likelihood selection, 24 candidate events are retained. The expected SM back- 
ground is 19.9±0.8(stat.)±2.9(syst.) events, computed using Monte Carlo simulations. 
Within the contribution from Z°Z°— >qqqq (10.4 events), more than 90% contain at least 
one Z° which decays into bb. The efficiency for mh = 95 GeV is (47.0±0.8(stat.)±1.6(syst.))%. 
The systematic errors for the signal selection efficiencies and the background estimates 
are assigned in the same way as described in 0]. The signal selection efficiency as a 
function of the h° mass is given in Table ^ which summarises the performance of all the 
SM channels. 

In order to verify the correct modelling of the b-tag algorithm in the high-multiplicity 
four-jet environment, the following cross-check has been performed. The measured charged 
particle tracks and calorimeter clusters of pairs of hadronic Z° decay events recorded in 
1998 at ^/s = rriz have been overlaid to form pseudo-events which have topologies very 
similar to e"'"e~— >Z°Z''— +qqqq events at LEP2 energies. In order to model the b tag prop- 
erly, the charged track parameters of the second Z° event have been adjusted such that 
the reconstructed primary vertices of both overlaid events match. The same procedure 
has been applied to simulated hadronic Z° decays. These pseudo-events have been passed 
through the same analysis chain as the high energy data and Monte Carlo. The same 
preselection as described above has been applied. The event b-tag likelihood Scvt is cal- 
culated using as inputs the two largest jet-wise b-tag variables - the same b-tag variables 
used in the h°Z° selection. The distribution for B^vt is shown in Figure ^ for overlaid data 
and overlaid Monte Carlo as well as their relative difference in the event tagging rate as 
a function of the cut on -Bcvt- The data and the Monte Carlo model agree to within 2%, 
independent of the cut on -Bevt- Hence, no additional systematic error has been assigned 
on the signal efficiency as a result of this cross-check. 

4.2 The Missing- Energy Channel 

The preselection of the missing-energy analysis is designed to enhance a signal charac- 
terised by two hadronic jets and missing energy in the presence of radiative Z° events, 
untagged two-photon events, and W"^W~— >qq£i/£ events. These Standard Model back- 
grounds also can produce final states with two hadronic jets and missing energy, but 
they can be rejected because their missing momentum is predominantly along the beam 



axis (in the case of the radiative T? events and untagged two-photon events), or there 
is a high-momentum lepton (in the case of W"'"W~— >qq£z/£ events). The preselection is 
described below. 

The initial requirements are intended to guarantee well measured events with a re- 
duced two-photon and Z/7* background: (1) the number of tracks satisfying the quality 
requirements used in |M must be greater than six and more than 20% of all tracks; there 
must be no significant energy in the forward detectors as described in no hits in the 
MIP plug detector with a significant charge deposition; the total transverse momentum 
P^s the visible mass m^\^ must satisfy 5 ■ P^jg + mvis > 100 GeV; the total visible 
energy E^i^ < 0.8 ■ ^i; and there must be less than 50% of the visible energy in the 
angular region | cos6'| > 0.90. The following requirements reduce the Z/7* contribution: 
(2) the polar angle of the missing momentum, ^miss, niust satisfy | cos^missl < 0.95 and 
the z-component of the visible momentum must satisfy |Pvisl < 35 GeV; (3) the tracks 
and clusters are grouped into two jets using the Durham algorithm, and the directions of 
both jets are required to satisfy | cos6'jet| < 0.95; (4) the acoplanarity angleQ of the two 
jets, (/>acop, must be larger than 5°. The following requirements reduce the contribution 
from four-fermion processes: (5) the missing mass mmiss must be consistent with mz'- 
60 GeV < mmiss < 120 GeV; and (6) no identified isolated leptons, as defined in 0, may 
appear in the event. 

The b-tag described in Section ^ is incorporated into the analysis by combining Bi and 
B2, the b-tagging discriminants of the two jets in the event, with four kinematic variables, 
I cos 6*111188 1, max(cos6'jet), "n^miss, and (/)acop, using the same likelihood technique as used at 
lower energies 0|. The high-pt lepton tagging has been removed from the calculation of 
the jet b-tagging discriminant variables Bi in order to avoid enhancing the W"*"W~— s>qq£z/£ 
background. Distributions of the kinematic and b-tag variables are shown in Figure |^ for 
OPAL data and SM background simulations. The signal likelihood is required to be 
larger than 0.60 for an event to be selected as a Higgs boson candidate. In Figure ^ 
the likelihood distribution and the mass distribution for the selected candidate events 
are shown for the data, SM backgrounds, and a simulated signal at m-h = 95 GeV. The 
reconstructed Higgs mass is evaluated using a kinematic fit constraining the recoil mass 
to the mass. 

The numbers of observed and expected events after each selection step are given^ in 
Table |1]. The selection efficiency estimated from the Monte Carlo for a 95 GeV Higgs is 
(35.4±0.9(stat.)±0.9(syst.))%. Ten events survive the selection, while 6.9±0.5(stat.)±0.6(syst.) 
events are expected from SM background processes. The systematic error evaluation is 
described in [Q. The detection efficiencies as a function of the Higgs boson mass are listed 
in Table |. 



4.3 The Tau Channels 

The tau channel selection consists of a preselection and tau lepton identification using an 
ANN, the details of which are described in The preselection requirements to select the 
signal with a low contamination of Z/7* events are: the event must be a high-multiplicity 

^ The acoplanarity angle, 0acop: of two vectors is 180° minus the opening angle between them in the 
plane transverse to the beam direction. 

In the calculation of the efficiencies and backgrounds in the missing-energy channel, a 2.5% relative 
reduction has been applied to the Monte Carlo estimates in order to account for accidental vetoes due to 
accelerator-related backgrounds in the forward detectors. 



multihadronic event P2|, |cos^^miss| < 0.95, the missing momentum pmiss niust be less than 



0.3 -a/s, the scalar sum of the transverse momenta of the particles in the event must exceed 
45 GeV, and at least one pair of oppositely-charged tau candidates must be identified. 

To select signal candidates, we use the two-tau likelihood of = p^pg_|_(f_!p^)(i_p^) , 

where Pj is the probability that the i^^ tau candidate originates from a real tau lepton. 
This probability is calculated from the shapes of the ANN output for signal and fake 
taus. The distribution of the ANN output for signal events was computed from Monte 
Carlo simulations. This analysis has been improved by using for the fake taus the distri- 
bution of the ANN output in hadronic Z° decay data collected in the calibration run at 
a/s ^91 GeV, which has a low fraction of events with real taus. This estimation of the 
fake tau ANN distribution reduces the systematic uncertainty on the fake tau rate. To 
pass the selection, an event must have Crr of at least 0.10. 

Next, the particles in the events are subdivided into two tau candidates and two jets. 
A 2C kinematic fit is applied using total energy and momentum conservation constraints, 
where the tau momentum directions are taken from their visible decay products while 
leaving their energies free. The x^-probability of the fit is required to be larger than 10~^. 
In events where both taus are classified as one-prong decays, the sum of the momenta of 
the charged particles assigned to the tau decays must be less than 80 GeV, in order to 
reduce backgrounds from Z^Z^^*^ ^fi^ fi^qq and Z^Z'^^*''— >-e"'"e~qq. 

The two final likelihoods, described in are then formed. One, /^(bbr+r"), is 
optimised for the h°Z°— >bbr"^r~ final state and makes use of the b-tag of Section 0. The 
other, /^(qqr+r"), is optimised for the h°Z°— »-r"'"r~qq process and does not use b-tagging. 
The following variables described in 0] are inputs to both likelihoods: -Rvis= E^is/^/s, 
|cos6'miss|, '(^TT, the logarithm of 7/34, the energy of the most energetic electron or muon 
identified in the event (if any), the angles between each tau candidate and the nearest 
jet, and the logarithm of the larger of the two 3C kinematic fit probabilities, in which the 
additional constraint comes from fixing either the tau pair invariant mass or the jet pair 
invariant mass to the Z° mass. We have introduced a new variable which takes advantage 
of the finite lifetime of the tau lepton, the impact parameter information of charged 



tracks belonging to the tau candidate, combined in a joint-probability calculation [33 



In addition to the above variables, the £(bbr+r~) likelihood uses the output of the b- 
tagging algorithm described in Section ^. An event is retained if £(bbr"'"r~) exceeds 0.92 
or £(qqr"'"r~) exceeds 0.88. Figure ^ shows the distributions of Crr, the joint impact- 
parameter probability, and the two likelihoods, £(qqr"^r^) and £(bbr"*"r~), for the data, 
simulated Standard Model backgrounds, and also for a simulated 95 GeV Higgs signal. 

The numbers of observed and expected events after each stage of the selection are given 
in Table 0, together with the detection efficiency for a 95 GeV SM Higgs boson, which 
is estimated to be (34.3 ± l.l(stat.)± 2.4(syst.))% after the final selection requirement. 
Three events survive the likelihood cut, to be compared to the expected background of 
4.0 ± 0.5(stat.) ± 0.9(syst.). The systematic errors are evaluated as in 

These results are confirmed by a separate analysis which uses a different technique to 
tag the tau candidates. Events are reconstructed as four jets using the Durham algorithm 
and tau candidates are sought in the four jets using a likelihood technique. The b- 
tagging algorithm of Section |^ is applied to the two remaining jets. This selection has a 
performance similar to that of the ANN method for tau tagging described above. The 
efficiency for a 95 GeV Higgs signal is evaluated to be (30.6 ± l.l(stat.) ± 1.6(syst.))%. 
Six events are retained, two of which are selected also by the analysis above, compared 



with the expected background of 3.0 ± 0.5(stat.)±0.6(syst.) events. 

4.4 The Electron and Muon Channels 

The preselection for this analysis is intended to enhance the qq topology, where i =e 
or /i. The requirements are listed below. 

1. The event must have at least six charged tracks, 1/34 > 10^"^, \P^[s\ < (-^vis — 0.5a/s) 
and E^is > 0.6^/s. 

2. At least one pair of oppositely charged leptons of the same flavour (e or /x) must be 
identified as described in 0. Figure ||(a) shows the distribution of the normalized 
dE/dx, an important ingredient of the electron identification. 

3. The events are reconstructed as two leptons and two jets. In the case of the muon 
channel, a 4C kinematic fit, requiring total energy and momentum conservation, is 
applied to improve the mass resolution of the muon pair and is required to yield a 

probability larger than 10~^. The invariant mass of the lepton pair is required to 
be larger than 40 GeV. The distribution of reconstructed dimuon masses is shown 
in Figure ||(b). 

Two likelihoods, one based on kinematic variables, /C, and one for b-tagging for the 
two jets of hadrons, B2jet, are calculated as described in [Q. The final signal likelihood is 
computed using /C and B2jet as inputs. The combined likelihood is required to exceed 0.2 
for the electron channel and 0.3 for the muon channel. Distributions of these likelihoods 
are shown in Figure |[ The signal selection efficiency for a 95 GeV SM Higgs boson is 
(55.0±0.9(stat.)d=l.l(syst.))% for the electron channel, and (64.9±0.9(stat.)±0.9(syst.))% 
for the muon channel. The numbers of observed and expected events after each stage of 
the selection are given in Table |l|, together with the detection efficiency for a 95 GeV SM 
Higgs boson. The selection retains three events in the electron channel and one in the 
muon channel. The total background expectation is 2.6±0.2(stat.)±0.5(syst.) events in 
the electron channel and 2.1±0.1(stat.)±0.4(syst.) events in the muon channel. 

4.5 Search for the MSSM Process h^Z^ with h^^A^A^ 

If 2mA < trih, the decay h*^— i>A'^A° is kinematically allowed and is the dominant de- 
cay of the h° in parts of MSSM parameter space. Dedicated, optimised searches for 
h°Z°— ^A'^A'^Z^ have not been performed. Instead, we evaluate the sensitivity of the h°Z° 
searches for processes with h°— i>A°A° by studying the efficiencies of the selections using a 
h'^Z'^— >-A°A'^Z°— i>bbbbZ° Monte Carlo. Non-bb decays of the A*^ are not considered. The 
efficiencies of the h°Z° selections for the process h°Z°— >-A°A°Z° with ttt-a = 20 GeV and 
mh = 70 GeV are summarised in Table |[ In general, the efficiency is found to increase 
with the h° mass. This procedure has the advantage of simplifying the statistical treat- 
ment and is sufficiently powerful to exclude regions in which the decay h'^^A^A'^^bbbb 
dominates. The expected backgrounds and the selected candidates are the same as for 
the results described above, since the same selections are applied. 
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Table 1: The h°Z" channels: the numbers of events after each cut for the data and the 
expected background, normalised to the data luminosity. The two-photon background, 
not shown separately, is included in the total background. The last column shows the 
detection efficiencies, forh°— > bb in the four-jet channel, forh°— > all in the missing-energy 
electron, and muon channels, and for Z'^h^— >T"'"T~(h°^ all) or Z'^h^— >qqT'^T~ in the tau 
channel, for a Higgs boson mass of 95 GeV with Standard Model branching fractions. 
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Table 2: Signal detection efEciencies for the searches for the SM Higgs boson, apphed 
to the processes with h''— >A°A° followed by A^—yhh. The efficiencies are quoted for 
mh = 70 GeV and = 20 GeV. The statistical errors due to the limited sizes of the 
Monte Carlo samples are 1-4%. The backgrounds and candidates for each channel are 
listed in Table ||. 



5 The A h search channels 

We search for the MSSM process of associated production e"'"e~— i>A°h° followed by the 
decays A°h°— >bbbb and A°h''^bbr"'"r~. If h°— >A°A'' is kinematically allowed, the process 
e"'"e~— >-h°A°— >A'^A°A° is also searched for in the bbbbbb final state. These selections differ 
from those for Z^h^ because the Z° mass constraint is no longer applicable. The selections 
therefore have been optimised to reject backgrounds differing kinematically from those 
important to the Z'^h° searches. 



5.1 The A^hO^bbbb Final State 

The same preselection requirements are applied as described in |Q. These are: 
1-3. The requirements 1-3 of the h^Z^ four-jet analysis of Section are used. 

4. The C-parameter must be larger than 0.45. 

5. Each of the four jets must contain at least six energy-flow objects and at least one 
charged track. 

6. The probability of a 4C flt, which requires energy and momentum conservation, 
must be larger than 10~^. 

For events passing the preselection, a likelihood technique is applied as described in [^. 
Seven input variables are used. Six of these variables are the same as those used in [^]: the 
four b-tagging discriminants Bi, one for each jet (see Section P), 7/34 in the Durham scheme 
and the event thrust. The seventh variable is (| cosGdijetl), the average absolute value of 
the cosines of the polar angles of the momenta of the two di-jet systems for the combination 
yielding the smallest difference in the two di-jet invariant masses after the 4C kinematic 
flt. This last variable replaces the mean |cos^jet| of the four jets used in 0] because 
it provides better discrimination against the W^W~ background. The distribution of 
(I cosGdijetl) is shown in Figure along with the flnal likelihood discriminant C^^. 

Candidate events are selected by requiring C^^ > 0.95. Table ^ shows the number of 
selected events together with the expectation from background processes and the signal 
selection efficiency for m,A= mh= 80 GeV, after each cut in the preselection and after the 



final cut on C^^. The detection efficiency for a Higgs signal with masses mp = mh= 80 GeV 
is estimated to be (48.4±0.7(stat.)±3.9(syst.))%. Eight candidate events are observed in 
the data, consistent with 8.0±0.5(stat.)±1.4(syst.) events expected from SM background 
processes. Six of the candidate events are common to those found in the SM four-jet 



channel of Section ^4.1| . The systematic uncertainties on the signal selection efficiencies 
and background estimates were determined using the same methods as described in 

Candidate Higgs masses are calculated from the measured jet momenta using the 4C 
ffi. Figure |TD|(a)-(c) shows the distribution of the sum of the reconstructed Higgs masses, 
Msuni = 'rrih'^ + 'rn^A, for all three possible di-jet pairings, separately according to the 
reconstructed Higgs mass difference, AM = \m^l^ — mjif^l. The resolution on the mass 
sum, Msum, is estimated to be approximately 3 GeV for M = 150 GeV. For mh = rriA, 
68% of the events have a reconstructed mass difference AM of less than 13 GeV. Since the 
four jets can be combined in three ways, and since the h° and A° cannot be distinguished, 
each candidate event enters at six points in the (mh,mA) plane. 

5.2 The AOhO^bbr+T- Final State 

The A°h°^bbr+r^ final state, where either the A° or the h° decays into a tau pair, has 
been searched for using the same technique as used in the SM tau channels described 
above. The final likelihood selection has been optimised for the MSSM process. 

The following variables are used as inputs to the calculation of the likelihood, £hA: -Rvis, 
|cos^jniss|, '^TT, the logarithm of 1/34, the energy of the most energetic identified electron 



or muon, the tau-track lifetime information as described in Section [4.3| , the average of the 
absolute values of the cosines of the polar angles of the reconstructed h° and A*^, which 
gives an estimate of the production angle of the two Higgs bosons, and the outputs of the 
b-tagging algorithm described in Section ^ for the two hadronic jets in the event. 

Candidate events are selected by requiring the final likelihood value to exceed 0.64. 
Table 1^ shows the number of selected events, the efficiency for m-h = rriA = 80 GeV and the 
background estimation, after each stage of the selection. The efficiency for a signal with 

= = 80 GeV is estimated to be (45.3±1.5(stat.)±2.3(syst.))%. Seven candidates 
are observed in the data, two of which are shared with the SM analysis of Section [4.3| . The 
number of events expected from SM background processes is 4.9 ± 0.6(stat.)±1.6(syst.). 
The distribution of the sum of the reconstructed Higgs masses, M^*"^ + Mj^'^, is shown 
in Figure |l3(d). Since the A° and h° cannot be distinguished, each selected event has 
two interpretations in the (mh, rriA) plane. Systematic uncertainties on backgrounds and 
efficiencies are evaluated as in Section |4.3| . 

The alternative jet-based analysis, described in Section [4.3| , is also applied in this 
search channel. This procedure gives an overall efficiency of (39.1±1.5)% for mh = itla = 
80 GeV, similar to that obtained above with the ANN tau identification, and the effi- 
ciency remains similar over the range of kinematically allowed values of mh and rriA- The 
choice of the ANN analysis has been made to optimise the expected limits. Five events, 
two of which are shared with the ANN analysis, meet the selection requirements of the 
alternative analysis, consistent with the expectation of 4.8 ± 0.6(stat.)±0.9(syst.) from 
SM background processes. 



5.3 The AOhO^AOA^A^^bbbbbb Final State 



Signal events in this decay mode are characterised by a large number of jets containing 
b-flavoured hadrons. Backgrounds are reduced through cuts on kinematic variables and 
by b-tagging. Signal events satisfy 2mA < and this analysis uses an ANN to optimise 
the sensitivity over the allowed mass region in the (mh,mA) plane. 

The analysis begins with a set of preselection requirements. Events must be selected as 



hadronic final states [^. The polar angle of the thrust axis, 9t, must satisfy | cos^tI < 
0.9. The charged tracks and clusters are grouped into six jets using the Durham ^9 
algorithm. Each jet is required to have at least one charged track and one electromagnetic 
cluster. Events must have at least 20 good quality charged tracks and at least 20 good 
quality electromagnetic calorimeter clusters, where the quality requirements are those 
used in [Q. To suppress the background from (Z/7)*, the value of yu is required to be 
larger than 0.0005 and the C-parameter to be larger than 0.0075. 

After the preselection, candidates are selected using two neural networks applied se- 
quentially, one combining kinematic and topological variables, and one for the b-tagging 
variables. The kinematic characteristics of signal events depend strongly on and rriA, 
and in order to maintain good sensitivity over a broad range of these masses, 14 variables 
with complementary discriminating roles are used as inputs for the first network. These 
are the event thrust, sphericity, oblateness, C and D parameters [^, the 2"^^ through 
the 6*^ normalised Fox- Wolfram moments, and ?/23, 1/34, I/45 and for the Durham 
scheme. Events were required to have a kinematic network output greater than 0.68. 

The final selection was made with an ANN optimised for b-tagging. If the A° is light, 
the two jets originating from the A° decay may be observed as a single jet; events then 
may seem to have only three jets. The events are therefore reconstructed both as three 
and six jets and a 15-input neural network has been trained using the b-tag variables 
described in Section ^, from the three and six jets in the event. Events are required to 
have a b-tagging network output greater than 0.92. The distributions of the two network 
outputs are shown in Figure |lT] for events passing the preselection. 

Five events pass all selection requirements (Table consistent with an expected back- 
ground of 8.7±1.0(stat.)±2.5(syst.) events. The signal efficiency for mA=20 and mh=70 
GeV is (45.4±2.2(stat.)±4.3(syst.))%. The systematic errors are evaluated similarly to 
those for the other channels 0] . The dominant contribution to the systematic error on the 
signal and background efficiencies comes from the uncertainty related to the b-tagging. 
The main components arise from uncertainties in the tracking performance of the detector 
and the b-hadron decay multiplicity |]35| . 



6 Limits 

6.1 Mass limit for the SM Higgs boson 

We summarise the efficiencies and numbers of expected SM Higgs signal and background 
events in Table |^, separately for 189 GeV and with the 183 GeV sample added. A 
total of 41 events at 189 GeV pass the selections while 35.4 ± l.l(stat.) ± 3.2(syst.) 
events are expected from the SM background processes. The reconstructed masses of the 
candidate events and the distribution expected for the combined SM backgrounds for the 
searches presented here, added to those from our ^/s ~ 183 GeV searches |Q, are shown 
in Figure |12[ A slight excess of events peaks at the Z" mass, which is consistent with a 
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Table 3: The A^h^ channels: effect of the cuts on data and the simulated background, 
normalised to the integrated luminosity of the data. The two-photon backgrounds, not 
shown separately, are included in the total background. The signal efficiencies are given 
in the last column for mii=mx=80 GeV in the h^A''— i>bbbb and the tau channels and for 
mh = 70 GeV and mA = 20 GeV in the h°AO^AOA°A°^bbbbbb channel 

statistical ffuctuation of the dominant Z^Z^ background. 

To derive a lower bound on the mass of the SM Higgs boson, we combine the results 
of the searches presented here with those of our earlier searches at ^/s 183 GeV [^. 
opal's earlier searches at energies of 172 GeV and lower have a negligible impact on the 
expected limit because of the limited kinematic reach of the lower-energy searches. The 
systematic errors on the background expectations and on the signal selection efficiencies 
have been treated following the method given in Reference [^. The confidence level for 
the signal hypothesis, CLg, is shown in Figure |13| (a). It has been computed using the 
weighted event-counting method described in Section 5 of Reference 0. 

The number of expected signal events and the upper limit on the production rate for 
signal events at the 95% confidence level (CL) are given in Figure |T3] (b) as functions of 
the Higgs boson mass hypothesis. A lower mass bound of 91.0 GeV is obtained at the 
95% CL, while the average expected limit from a large number of fictitious experiments 
assuming the background-only (zero signal) hypothesis is 94.9 GeV. The probability to 
obtain a limit of 91.0 GeV or less in an ensemble of background-only experiments is 4%. 
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Table 4: Detection efEciencies in percent and numbers of expected SM Higgs boson events 
(in parentheses) , at y/s = 189 GeV, for each search channel separately, as a function of the 
Higgs boson mass. The last two rows show the sum of expected SM backgrounds, the error 
on the background sum, and the numbers of candidates, for each channel separately. The 
last two columns show the total numbers of expected events in all channels for the present 
search at y/s = 189 GeV, and also summed with the results obtained at 183 GeV 



6.2 Limits in the MSSM parameter space 

The searches presented in this pubhcation at ^/s ~ 189 GeV are combined with previous 
OPAL Higgs searches 0, ||, |, ^ at ^/s between mz and 183 GeV. We present 

95% CL hmits in the MSSM parameter space for the constrained MSSM with six param- 
eters, mo, M2, A, /i, tan/3 = V2/V1 and ttt-a in addition to those of the SM. The definitions 
of these parameters can be found in The consistency of the MSSM Higgs predictions 
with our observed data has been calculated for a large sample of models in the space 
spanned by these six parameters. The scanning strategy is described in detail in |^. 

Two "benchmark" situations (Scan A of 0) are first considered, corresponding to no 
mixing and maximal mixing in the scalar-top sector (A = and a/6 TeV, respectively), 
with 5 < rriA < 2000 GeV and 0.7 < tan (3 < 50. The remaining free parameters are fixed 
to "universal" values: the soft SUSY breaking mass parameters at the electroweak scale 
{e.g., mq, the "left-up" scalar quark mass) are set to 1 TeV, M2 is setQ to 1.63 TeV and /i 
is fixed to —100 GeV. The top quark mass is fixed at 175 GeV. The exclusions are shown 
in Figure [14| in the (mh,mA) plane, the (mh,tan/?) plane|, and the (rriA, tan/3) plane. Both 



cases of scalar top mixing are included. If for either choice of the mixing a parameter set 
is not excluded, it is shown as unexcluded. 

For tan/? > 1, lower mass limits of mh > 74.8 GeV and > 76.5 GeV are obtained, 
while the expected limits are mh > 76.4 GeV and mA > 78.2 GeV. In the case of no scalar 
top mixing, we exclude the range 0.72 < tan/3 < 2.19, while no tan/3 is excluded in the 
case of maximal scalar top mixing. For smaller mtop the excluded region of tan f3 becomes 
larger, while for larger mtop, it becomes smaller. If tan/3 is allowed to vary between 0.7 

^ In References Q and |^ it was mistakenly stated that M2 was set to 1.0 TeV for Scan A. The actual 
value of M2 used was 1.63 TeV in those papers. 

^The corresponding Figure 19(c) of Reference B mistakenly showed the 99% CL exclusion contour. 



and 1.0, the range of possiblilities in the model is larger and a small, unexcluded region 
appears, visible in Figures |1^ (b), (c) and (d), with mh~70 GeV, and mA<10 GeV. 
This region corresponds to models for which the dominant decay of the h° is into A^A*^, 
while A° decays into bb are kinematically impossible, rendering the analyses of this paper 
inefficient. Nearby regions are excluded because the decay h"— i>bb proceeds, although 
with a reduced branching fraction. 

In a more general scan (Scan C of [P]), all six parameters are varied independently 
within ranges motivated by theory, with the top quark mass taken at 165, 175 and 
185 GeV, which allows for about ±2a of the current measurement error. The ranges for the 
parameters of this scan are < mo < 1 TeV, < M2 < 2 TeV, —2.5 ■ mo < A < 2.5 ■ mo, 
— 1 < /X < 1 TeV, 5 < rriA < 2000 GeV and 0.7 < tan/5 < 50. For some parameter 
sets, the heavier CP-even Higgs (H°) is light enough to be produced in the H°Z° process, 
and in those cases the h°Z° searches are considered sensitive to it, thereby extending the 
excluded area. Parameter sets giving rise to chargino or neutralino masses |^ or scalar 



top masses already excluded by OPAL searches are considered excluded here. Models 
which give rise to large Z°— i>h°Z* or Z*^— >-h°A° cross-sections incompatible with the mea- 
sured Z° decay width (see [0) are also considered excluded. The results for this scan are 
shown in Figure |15|. The dark area is excluded at the 95% CL. The unexcluded region for 
mh~70 GeV and mA<10 GeV has become larger in this more general scan, but it is still 
limited to the region of low tan f3. 

For large values of A and /i, it is possible for the MSSM Lagrangian to have minima 
which have non-zero vacuum expectation values for the t fields which break charge and 
colour symmetry . One may identify parameter sets for which the Lagrangian has such 



a minimum 



J3[, but the criteria for excluding such sets because they are inconsistent with 
observation are substantially modified if it is possible for the electroweak minimum to be 
only a local minimum and the tunneling rate to the charge- and colour-breaking (CCB) 
minimum is longer than the age of the universe. A specific model including the effect of 
tunneling is available out of a number of different possibilities . A simple, approximate 



criterion to avoid CCB minima is 42 



A^ + 3fi^ < x{ml + ml^] 



where m^^^ and mi^ denote the left- and right-handed scalar top masses and x ~3. For the 
specific calculation that includes the possibilities of false vacua which tunnel only slowly 
to CCB vacua, this bound is shown to be modified to x ~ 7.5. We show in Figure |TB| the 
regions in m^, mAand tan/? in our general scan which we do not exclude but which lead 
to Lagrangians with CCB minima according to the looser criterion with x = 7.5. With 
this CCB criterion applied, absolute mass limits mA>76.0 GeV and mh>72.2 GeV are 
derived for tan/5>l at the 95% CL. 



7 Conclusions 

A search for neutral Higgs bosons has been performed based on the data collected at 
^/s ~ 189 GeV with an integrated luminosity of approximately 170 pb~^. Searches have 
been performed for the Standard Model process e"'"e~— >H°Z° and the MSSM processes 
e"'"e~^h''Z'', A°h°. The search channels are designed to detect bb and r+r^ decays of 
the Higgs bosons, and the process h°— >-A°A° is also considered. No significant excess of 
candidates is observed in the data beyond the expected Standard Model backgrounds, and 



we derive the following limits at the 95% confidence level. For the SM Higgs boson, we 
obtain a lower mass bound of 91.0 GeV. In the MSSM, we obtain the limits mh> 74.8 GeV 
and mA> 76.5 GeV assuming tan/3>l, the mixing in the scalar top sector to be either zero 
or maximal, and the soft SUSY-breaking masses are 1 TeV. In the general scan of MSSM 
parameters, excluding parameter sets which result in Lagrangians with charge- and colour- 
breaking minima, we derive absolute mass limits of mh > 72.2 GeV and > 76.0 GeV 
for tan/3> 1 at the 95% confidence level. 
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Figure 1: B-tagging performance and modelling, (a): the distribution of the output of 
the b-tagging algorithm, B, for jets in calibration data taken at ^/s = mzo, compared to 
the Monte Carlo expectation. The data distribution is given by the points, with error 
bars smaller than the plot symbols. The open histogram shows the distribution of B for 
b-flavoured jets, and the dark (light) grey histogram shows the contribution from c (uds) 
flavoured jets, expected in a Monte Carlo simulation. Inset: The b-tagging performance 
in the present work for hadronic jets in Z° decay compared with our previous version 
(b): the relative difference of the tagging rates for jets opposite to b-tagged jets in Z° 
calibration data and Monte Carlo as a function of the likelihood cut. The points are 
statistically correlated because the tagging rate accumulates events on one side of the 
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Figure 2: Input variables for the four-jet channel likelihood selection: (a) the fit proba- 
bility of a 6C kinematic fit which requires energy and momentum conservation and both 
dijet masses to be equal to m-^; (b) the logarithm of the jet resolution parameter 7/34; 

(c) the difference between the energies of the jets with the highest and lowest energies; 

(d) the second largest jet b-tag. OPAL data are indicated by points with error bars, 
four-fermion backgrounds by the light grey histograms, and two-fermion backgrounds by 
the dark grey histograms. The estimated contribution from a 95 GeV Higgs is shown with 
dashed histograms; it has been scaled by a factor of 50. 
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Figure 3: Four-jet channel: the hkehhood distribution (a) and the reconstructed candidate 
masses (b). OPAL data are indicated by points with error bars, four-fermion backgrounds 
by the light grey histograms, and two-fermion backgrounds by the dark grey histograms. 
Also shown is the contribution expected from a 95 GeV Higgs boson (hatched histograms). 
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Figure 4: Four-jet channel: (a) the Bevt distribution for overlaid pairs of hadronic Z° 
decay events (described fully in the text). OPAL data are indicated by points (the error 
bars are smaller than the points), and the Monte Carlo simulation by the histograms with 
the contributions from final states containing four b-quarks (light), two b-quarks (dark) 
and only light flavoured quarks (open). The discontinuous behaviour in the central bins 
is due to the discrete binning of the two b-tag variables used as inputs to the likelihood, 
(b) The relative difference between OPAL data and the Monte Carlo simulation for the 
event b-tagging rate as a function of the cut on B^vt for pairs of overlaid Z° decays. 
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Figure 5: Input variables to the missing-energy channel likelihood selection: (a) the 
acoplanarity angle; (b) the invariant mass recoiling against the hadronic system; (c) the 
cosine of the polar angle of the missing momentum; (d) the larger cosine of the two jet 
polar angles; (e) and (f) the b-tags of the two jets. OPAL data are indicated by points 
with error bars, four-fermion backgrounds by the light grey histograms, and two-fermion 
backgrounds by the dark grey histograms. The estimated contribution from a 95 GeV 
Higgs boson is shown with dashed histograms; it has been scaled up by a factor of 20. 
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Figure 6: Missing-energy channel: (a) the hkehhood distribution and (b) the reconstructed 
candidate masses. OPAL data are indicated by points with error bars, four-ferniion 
backgrounds by the hght grey histograms, and two-fermion backgrounds by the dark grey 
histograms. Also shown is the contribution expected from a 95 GeV Higgs boson (hatched 
histograms). 
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Figure 7: Tau channels: (a) distribution of the two-tau hkehhood as input to the hkeh- 
hood selection; (b) distribution of the joint impact parameter significance as input to the 
likelihood selection; (c) distribution of the qqrr likelihood; (d) distribution of the bbrr 
likelihood; OPAL data are indicated by points with error bars, four-fermion backgrounds 
by the light grey histograms, and two-fermion backgrounds by the dark grey histograms. 
Also shown is the contribution expected from a 95 GeV Higgs boson (dashed histograms 
in (a) and (b), hatched histograms in (c) and (d)). In (a) and (b) the signal is scaled by 
a factor of 100. 
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Figure 8: Electron and Muon channels: (a) The distribution of the normalised dE/dx, the 
difference between the ionisation energy loss of a track and that expected for an electron, 
divided by the measurement error, for tracks in the high-energy data, (b) Reconstructed 
invariant mass of identified muon pairs, where the signal is expected to contribute near 
TTT'/ifi = tTLio. (c) The electron channel likelihood distribution, and (d) the muon channel 
likelihood distribution. OPAL data are indicated by points with error bars, four-fermion 
backgrounds by the light grey histograms, and two-fermion backgrounds by the darker 
grey histograms. The contribution from a 95 GeV Higgs boson signal is shown with 
dashed histograms in (a) and (b) and with hatched histograms in (c) and (d). It is scaled 
by a factor of 100 in (a) and (b). 



OPAL 



!/5 



in 
o 



225 
200 
175 
150 
125 
100 
75 
50 
25 


10 



OPAL data 
□ 2-fermion 
I I 4-fermion 
--- hA-signal x 100 
(mA= iTij, = 80 GeV) 



(a) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

< I COsQdijet > 



OPAL data 
I I 2-fermion 
I I 4-fermion 

hA-signal (m^= m^ = 80 GeV) 



(b) 
cut 




0.8 1 

Ah-Likelihood 



Figure 9: (a) The distribution of (| cos6dijct|) for the A°h°^bbbb channel, (b) The dis- 
tribution of the selection likelihood in the A°h°— i>bbbb channel. OPAL data are indicated 
by points with error bars, four-fermion backgrounds by the light grey histograms, and 
two-fermion backgrounds by the dark grey histograms. Also shown is the contribution 
expected from a Higgs boson signal with = mp^ = 80 GeV (other MSSM parameters 
are set to their "benchmark" values with maximal scalar top mixing, as described in Sec- 



tion p^ ) as dashed (a) and hatched (b) histograms. In (a) the signal is scaled by a factor 
of 100. 
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Figure 10: The sum of the reconstructed Higgs masses, m^^^ + mJi^^, for (a) the Ah—>-4b 
channel with the di-jet pairing combination which yields the smallest mass difference, 
AM = |m^'^ — mj^'^'^l, (b) the Ah—^4b channel with the medium AM combination, (c) the 
Ah^4b channel with the maximum AM combination, and (d) the Ah^ r+r~bb, bbr+r~ 
channel. OPAL data are indicated by points with error bars, four-fermion backgrounds by 
the light grey histograms, and two-fermion backgrounds by the darker grey histograms. 
Also shown as dashed histograms are the contributions expected from a Higgs boson signal 
with mh = rriA = 80 GeV (other MSSM parameters are set to their ^'benchmark" values 
with maximal scalar top mixing, as described in Section |6'.2| )- 
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Figure 11: e+e"->A°h°->A°A°A°^bbbbbb channel: The distribution of the ANNs with 
kinematic variables (a) and with b-tagging variables (b), as described in the text. The 
light grey histograms indicate Standard Model qq background, and the dark grey his- 
tograms show the distributions of Standard Model four-fermion background. The dashed 
histograms show the distribution from a signal Monte Carlo with Higgs boson masses of 
m}^=70 GeV and mp.^=20 GeV (other MSSM parameters are set to their benchmark" 
values with maximal scalar top mixing, as described in Section scaled by a factor of 
10. 
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Figure 12: The mass distribution for the selected events in all h^Z" channels combined 
(points with error bars), and the expected background. The data and expected back- 
grounds have been combined for the =183 GeV and = 189 GeV samples. The 
expected mass distribution assuming the production of the SM Higgs boson with a mass 
of 91 GeV (95 GeV) is added on top of the background and shown with a dashed (dotted) 
histogram. 
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Figure 13: Standard Model Z'^U.^ limits: (a) The conGdence level for the signal hypothesis 
as observed (solid line) and expected on average for background-only experiments (dashed 
line), as functions of the SM Higgs boson mass, (b) The limit on the production rate for 
the SM Higgs boson at 95% CL (solid line) and the number of expected signal events 
(dashed line) as functions of the Higgs boson mass. The 95% CL lower mass limit on 
the SM Higgs bosons is set at the point where the solid and dashed curve intersect: 
mn > 91.0 GeV. 
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Figure 14: The MSSM exclusion for the "benchmark 
tion \6.2^ . Excluded regions are shown for (a) the (mh 
(rrih, ^a) plane for tan/? > 0.7, (c) the (mh, tan/?) plane, and (d) the {rriA, tan/?) plane. 
The black area is excluded at the 95% CL. The grey areas in (a), (b) and (c) are theoret- 
ically inaccessible. In (b), because of the wider permitted range of tan f3, more {m^, mp^ 
points become available, and a small, unexcluded region appears at mh~70 GeV, with 
mA<10 GeV; it can also be seen in (c) and (d). Shown in (d) is the excluded region for 
no scalar top mixing. In all figures, the black region is excluded for all values of the scalar 
top mixing. In (d), the grey area is excluded for the case of no scalar top mixing. The 
dashed lines indicate the boundary of the region expected to be excluded at the 95% CL 
if only SM background processes are present. The dotted line in (d) is the expected limit 
for the case with no scalar top mixing. 
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Figure 15: The MSSM exclusion for the general scan described in the text of Section \B.2i . 
Excluded regions are shown for (a) the (mh, mp^) plane for tan/? > 1, (b) the (mh, mj^) 
plane for tanP > 0.7, (c) the (mh, tan/3) plane, and (d) the {mp^, tan/3) plane. All exclusion 
limits are at the 95% CL. The black excluded without applying a CCB criterion 

(described in the text). The grey hatched excluded when the CCB criterion is 

applied with x = 7.5. The grey areas in (a), (b) and (c) are theoretically inaccessible. 



